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Abstract 

m : 

We study the constraints on flavour violating terms in low energy SUSY coming 
from several processes as li — > Ijj, li — > Ijljlj and /i — > e m Nuclei. We show that 
{Sj ' a combined analysis of these processes allows us to extract additional information 

with respect to an individual analysis of all the processes. In particular, it makes 
possible to put bounds on sectors previously unconstrained by k — > Ij-y. We perform 
the analysis both in the mass eigenstate and in the mass insertion approximations 
clarifying the limit of applicability of these approximations. 
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1 Introduction 

M-l, 

Neutrino oscillation experiments have established the existence of lepton 
family number violation processes. So, as a natural consequence of neutrino 
oscillations, one would expect flavour mixing in the charged lepton sector. This 
mixing can be manifested in rare decay processes such as /i — > cy, r — > /17, 
etc. However, if only the lepton yukawa couplings carry this information on 
flavour mixing, as in the Standard Model with massive neutrinos, the expected 
rates of these processes are extremely tiny being proportional to the ratio of 
masses of neutrinos over the masses of the W bosons [1] . These values are very 
far from the present and upcoming experimental upper bounds that we can 
read in Table 1. 

In a supersymmetric (SUSY) framework the situation is completely different 
[2]. For instance, the supersymmetric extension of the see-saw model [3] pro- 
vides new direct sources of flavour violation, namely the possible presence of 
off-diagonal soft terms in the slepton mass matrices and in the trilinear cou- 
plings [4]. In practice, flavour violation would originate from any misalignment 
between fermion and sfermion mass eigenstates. 

One of the major problems of low energy supersymmetry is to understand 
why all LFV processes are so suppressed. This suppression imposes very se- 
vere constraints on the pattern of the sfermion mass matrices which must be 
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Present and Upcoming experimental limits on various leptonic processes 
[5,6,7,8,9,10] 

either very close to the unit matrix in the flavour space (flavour universality) or 
almost proportional to the corresponding fermion mass matrices (alignment). 
Both universality and alignment can be either present as a kind of "initial" 
conditions or as a result from some dynamics of the theory. 
Given a specific SUSY model, it is possible, in that context, to make a full 
computation of all the FCNC (and possibly also CP violating) phenomena. 
This is the case, for instance, of the constrained minimal supersymmetric stan- 
dard model (CMSSM) where these detailed computations led to the result of 
utmost importance that this model succeeds to pass unscathed all the se- 
vere FCNC and CP tests. However, given the variety of options that exists in 
extending the MSSM (for instance embedding it in some more fundamental 
theory at larger scale), it is important to have a way to extract from the whole 
host of FCNC and CP phenomena a set of upper limits on quantities which 
can be readily computed in any chosen SUSY frame. Namely, one needs some 
kind of model-independent parameterization of the FCNC and CP quantities 
in SUSY to have an accessible immediate test of variants of the MSSM. 
The best parameterization of this kind that has been proposed is in the frame- 
work of the so-called mass insertion approximation (MI) [11]. One chooses a 
basis for the fermion and sfermion states where all the couplings of these par- 
ticles with neutralinos or with gluinos are flavour diagonal, while the flavour 
changing is exhibited by the non-diagonality of the sfermion propagators. De- 
noting with Ajj the off-diagonal terms (in flavour space) in the sfermion mass 
matrices, the sfermion propagators can be expanded as a series in terms of 
5ij = Aij/rh 2 where m is an average sfermion mass squared. As long as 5^ 
is much smaller than one, we can just take the first term of this expansion 
and then the experimental information concerning FCNC and CP violating 
phenomena translates into upper bounds on these 5ij. 

Obviously, the above mass insertion method is advantageous because one does 
not need the full diagonalization of the sfermion mass matrices to perform a 
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test of the considered susy model in the FCNC sector. 

Pioneering studies [12,13] considered only the photinos and gluinos mediated 
FCNC contributions because a complete inclusion of the neutralino and chargino 
sector contributions would require a complete specification of the model. On 
the contrary, the spirit of their study was to provide a model-independent way 
to make a first check of the FCNC impact on classes of SUSY theories. 
Subsequently, general studies in the leptonic sector [14] included the chargino 
and neutralino contributions using a generalization of the MI, that we call GMI 
[15], which consists in an approximation where the gaugino-higgsino mixings 
are also treated as insertions in the propagators of the charginos and neutrali- 
nos inside the loop. With this additional insertion approximation, it is not 
necessary to fix a particular scenario to analyze the dependence on the many 
mass parameters. However, in the lepton sector, interferences among ampli- 
tudes are generically present in the models and they would affect the limits on 
the 5ij. This happens, for instance, for 5 RR due to a destructive interference 
between the bino and bino-higgsino amplitudes. Moreover r — > /ry and r — > &-/ 
experimental bounds are not so stringent as in fi — > e'j so that we can have 
not so tiny insertions independently from cancellations. 

In this context, the target of this work is to study all the possible constraints 
on flavour violating terms in low energy SUSY coming from several processes 
as U — > Zj7, hi — > Ijljlj and ji — * e in Nuclei and, subsequently, to clarify 
which is the limit of applicability of the MI and of its generalization. 
The approach is to fix a specific model of known spectrum (CMSSM) and 
to compare, systematically, the predictions of the full computation [16] with 
respect to the approximate, MI and GMI, computations [18,14]. 
In particular, we focused on the RR sector where strong cancellations make 
the sector unconstrained. We show that these cancellations prevent us from 
getting a bound in the RR sector both at the present and also in the future 
when the experimental sensitivity on Br{li — > Ijj) will be increased. So, being 
our aim to find constraints in the RR sector, we examined other LFV pro- 
cesses as U — > Ijljlj and /i — > e in Nuclei. 

Finally we took into consideration the bounds on the various double mass in- 
sertions 5£ 3 B 5g D , with A, B,C,D = L, R. These limits are important in order 
to get the largest amount of information on SUSY flavour symmetry breaking. 



2 LEPTON FLAVOUR VIOLATION IN SUPERSYMMETRY 

The study of lepton flavour violation in SUSY scenarios is one of the most 
promising subjects in low energy supersymmetric phenomenology, in that it 
shows, if observed, a clear signal of physics beyond the Standard Model [19,20]. 
The source of LFV is the soft supersymmetry breaking lagrangian which, 
in general, contains a too large number of FV couplings, so that the pre- 
dicted branching ratios for LFV processes usually exceed phenomenological 
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constraints [21] : this is a typical example of the SUSY flavour-problem. 
The usual way to solve this problem is to consider Lagrangians which result 
from models that break SUSY in a flavour blind manner, as in mSUGRA or 
Anomaly mediated supersymmetry breaking (AMSB), etc [22]. 
Even in this case, in general, flavour is violated in the Lagrangian at the weak 
scale. For instance, LFV can be induced by the existence of new particles at 
high scale with flavour violating couplings to the SM leptons (as right handed 
neutrinos in a see-saw model [12]) or the presence of new Yukawa interactions, 
as in superGUTs theories [23,24]. In this last case the flavour violation is com- 
municated to low energy fields through renormalization effects. 
On the other hand, in models based on supergravity or superstring theories, 
nonuniversal soft terms are generically present in the high scale effective 
Lagrangian [25]. 

In models with flavour symmetry imposed by a Froggatt-Nielsen mechanism, 
flavour violating corrections to the soft potential could be large [26,27,28]. 
Irrespective to the source of these FV entries, our approach is to assume LFV 
at low energy and to try to bound the FV terms present in slepton mass ma- 
trices. 

The processes we are going to study are Zj — > Ijj, Zj — > Ijljlj and /i — > 
e in Nuclei that are mediated, at one loop level, by neutralinos, charginos and 
sleptons. The relevant interaction Lagrangian for the considered processes is: 

C = k (C t R AX P R + C^ AX P L ) Xa + h (N« X P R + N^ AX P L ) x° A lx. (1) 

We can always work in the basis where the charged lepton masses and 
the gauge couplings are flavour diagonal. In this basis, in general, the slepton 
mass matrix is not diagonal and its off-diagonal entries induce the LFV. 

Ft L R L 

Our aim is to use the couplings C iA ' x and N iAX , which are functions of the 
FV entries, to constrain the structure of the slepton mass matrix. 



3 MASS INSERTION APPROXIMATION 



In the spirit of the mass insertion approximation, we treat the off-diagonal 
elements of the slepton mass matrix as insertions. Our convention for the 
slepton mass matrices is: 



PA) 



m 2 L (l + 5 LL ) (A- fi tan (3)mi + m L m R 5 LR 

(A - n tan (3)mi + m L m R S LR m R (l + 5 RR ) 



) 









where mi and m R are, respectively, the average masses of the L and R sleptons 
and A = am with a ~ 0(1). 

The MI now corresponds to a development of the slepton propagators around 
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the diagonal with the average slepton masses, m 2 L and m R . In practice, we are 
working in the basis of leptons, neutralinos and charginos mass eigenstates 
and slepton weak eigenstates. 

The FV is parametrized through a flavour violating mass insertion in the 
virtual slepton line. In this basis, the gaugino couplings are flavour- diagonal 
and so the Lagrangian in eq.l takes the form: 

-C = iI¥a(n^P r + N^P L ) k + IWa{K^Pr + KPPl) k 

+ f$n(ct$PR + C^Pl) h + h.c, i = e,iM,r (2) 

where the coefficient C^q and Ng® (with B,C — (L, R)) are: 

-92{0 R ) A1 , C LR _-__(0 L)il2 
Nit ] = -^(On)a2 - ^(OnU , N$ = V2 9l (0 N ) Al 

where Ol,r and On are hermitian matrices that diagonalize the chargino and 
neutralino mass matrices, respectively. 

3.1 Neutralino-chargino sector: GMI approximation. 

In the slepton mass matrix, we were able to factorize the source of LFV, 
namely the mass insertions 5^. 

We would like to have a similar tool in the chargino-neutralino sector too. In 
fact, it is difficult to understand the dependence of physical quantities on the 
elements of the chargino-neutralino mass matrices. For instance, in the GMI 
approximation we treat both the off-diagonal elements (flavour violating or 
not) of the slepton mass matrices and the off-diagonal terms (flavour conserv- 
ing) of the chargino and neutralino mass matrices as mass insertions. 
In order to understand the validity of this other kind of approximation let us 
consider the chargino mass matrix M c : 



(W R H 2R ) 



M 2 



y \^2mw sin (3 fi 







I 


'we) 







+ h.c, 



diagonalized by 2 x 2 hermitian matrices Ol and R as O r MqO\ = (Mc)diag- 
Now, assuming that M 2jl , //, |M 2jl ±/i| ~ A » Mz, we obtain as approximate 
mass eigenstates W~ and H~ with masses M 2 and //, respectively. Actually, 
corrections to the spectrum start from the order 0(Mf/A 2 ), while the off 
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diagonal elements of the Or,Ol matrices are of the form Mw\M2Sp(cp) + 
licp(sp)\/\M% — fi 2 ] where sp = sin(3 and eg = cos j3. The approximations made 
follow naturally from scenarios like m-SUGRA in which a large fi term is re- 
quired in order to get correct electroweak symmetry breaking (see for instance 
table II in [17]). The neutralino spectrum can be analyzed in a very similar 
way. It turns out that a "bino-like" LSP can very easily have the right cosmo- 
logical abundance to make a good dark matter candidate, so from this point 
of view the large fi limit may be preferred. 

In this approximation the amplitudes of the analyzed processes have two kinds 
of contributions: one without off-diagonal MI in the gauginos propagator, i.e. 
pure B and W, the other with one MI mixing B — H, W — H, proportional 
to Mw- The approximation consists in keeping at most one flip insertion, 
only the terms of 0{M w /M i ,M w / fi), so we are neglecting terms of order 
0(M^ v /m 2 usy ). In this basis both charginos, neutralinos and sleptons are in 
the flavour-eigenstates and the expressions of the amplitudes are very easy to 
understand. 



4 mSUGRA spectrum 

The aim of our analyses is to bound all the leptonic Si/s and to establish 
the limit of applicability of the approximate methods presented above. 
In this context, we prefer to work in a well defined model (m-SUGRA) to re- 
duce the number of free parameters and to simplify the physical interpretation 
of the results. 

Let us recall the constraints arising in mSUGRA, where the universality as- 
sumption reduces the parameters at the Planck scale to a common scalar 
mass, m , a common gaugino mass, Mi/ 2 , and a universal A = am term with 
a ~ 0(1). 

At the low energy scale, after the RGE running, the parameters Mi, Mi, mi 
and rriR are obtained as follows: 

M^^^Im^Mu) 
Ui{Mu) 

m 2 L (m w )-m 2 L (M u ) + 0.5M 2 (Mu) + 0.04M 1 2 (M [/ ) 

m 2 R (m w )^m 2 R (M u ) + 0.15M 2 {Mu) 

where Mjj is the unification scale and the mSUGRA constraints are satisfied 
when Mi(Mjj) = M 2 (M C7 ) = M 1/2 and m R {Mu) = m 2 L (Mjj) = m\. 
A very important constraint in mSUGRA comes from the radiative elec- 
troweak breaking condition that requires a fine-tuned fi in order to fulfill the 
minimum condition: 
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. 2| Mi 2 1 + 0.5 tan 2 f3 lrf 0.5 + 3.5 tan 2 j3 
I" I + -f = + t.a,^-l ' 

so, the spectrum of the model is completely known given the following set of 
parameters: m , Mi/ 2 , a, tan/3, 



5 LFV processes: Z; — > Zj7 , /; — > 3/,- , /i — > e m Ti. 

Let us discuss the branching ratios of the LFV rare processes as k — > Zj-y. 
The amplitudes of the processes take a form 

T = m k e x Uj{p - q)[iq u <J Xv (A L P L + A R P R )\ui(p) 

where p and q are momenta of the leptons 4 and of the photon respectively. 
The chirality flip of this transition is the reason of the appearance of the 
factor in the operator. The branching ratio of U — > Ij'y is given by 



BRjlj^l^) _ 487r 3 a„ ij 2 {j 2 a 3 



A1\ 2 + \A^) ~ tan/3 
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with m a typical susy scale. In SUSY, this chirality flip can be implemented 
in the external fermion line or at Yukawa vertex or in the internal sfermion 
line through a flavour conserving FC mass insertion Aff L = (A — /xtan/3)mj. 
Each coefficient in the above formula can be written as a sum of two terms, 

Al,r = A n LR + Al R 

where A^ R and A C L R stand for the contributions from the neutralino loops 
and from the chargino loops respectively. 

Finally, we consider the — > 3lj and fi ^ e in Ti processes. Both of them 
get contributions from penguin-type diagrams (with photon or Z boson ex- 
changes) and from box-type diagrams. However the 7 penguin-type contribu- 
tion is enhanced by a tan j3 factor with respect to the other contributions so 
it dominates and one can find the simple theoretical relations: 

Br(h -> Ijljlj) ~ 7 x lO' 3 BR(k -> Z j7 ), 

£r(// -> einTi) ~ 6 x 10~ 3 Si?(// -> ey). 

Now, imposing the actual experimental upper limits on the above LFV pro- 
cesses, we get the upper bounds on the S^s from each process. In particular, 
the ratios among the S^s upper bounds from different processes are: 

cfi^e in Ti a r^ey 

°21 — 4 °21 
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We note that U — > Zj'y give the more stringent bounds on the 5i/s. However, 
these relations are no longer true in special regions, where, strong cancella- 
tions reduce Br(li — > ^7) by several order of magnitude and destroy the above 
relations (see the discussion about 5 RR bounds in the next section). 
In this paper, we will consider only the contributions from neutralino and 
chargino sectors. However, Higgs bosons (h°, H°, A ) are also sensitive to 
flavour violation and mediate processes such as \x — > e in Nuclei [29], r — > 3/i 
[30] or t — > /it] [31] or r — > fi(e)j [33]. The amplitudes of these processes 
are sensitive to a higher degree in tan/5 than the chargino/neutralino ones 
(the BRs grow as (tan/3) 6 , though they are suppressed by additional Yukawa 
couplings) and thus could lead to large branching fractions at large tan (3 [32]. 



6 Bounds on 5 lj from k — > ^7 , k — > 3lj , /x — > e m Nuclei 

The approach we follow to put bounds on the various is to consider 
only one mass insertion contributing at a time for the -Br(Zj — > ^7). Each off- 
diagonal FV entry in the slepton mass matrices is put equal to zero, except 
for the FV insertion we are interested to constrain. 

Now, what we mean by full computation is to diagonalize the slepton mass 
matrix (with only one 5^ term) and to use the experimental Br(li — > ^7) 
limits to impose constraints on each insertion type. 

Imposing that the contribution of each S lj does not exceed (in absolute value) 
the experimental bounds, we obtain the limits on the 5 u 's, barring accidental 
cancellations. 1 . 

Hence forward, to be as clear as possible, we will speak in the MI language. 

6.1 Bounds on 5ll 

The amplitudes proportional to 8ll receive both U(l) and SU(2) type 
contributions. In the first case we can have a pure B exchange, with chirality- 
flip implemented in the external (internal) fermion (sfermion) line or at yukawa 
vertex from B — H°. In the second case we have W and W — H exchange both 
for charginos and for neutralinos. However in the W case, because the SU (2) 
Gauge fields don't couple to right-handed fields, we can't make a chirality 
flip in the internal sfermion line. When the chirality is flipped in the external 
lepton line, the amplitude does not contain the \i mass term and is not tan (3 
enhanced. 

In the following we give the expression for the amplitude A^ 2 = (Al~ 2 )su(2) + 
(A^ 2 )u(i) °f the li — > Ijj processes in the GMI approximation: 



In SUSY GUT theories there are possible correlations between the hadronic and 
leptonic FV effects [34,35]. 
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f A ij\ = ^2 A « 
V W 577(2) 47]- LL 



ml (M|-/i 2 ) mi 



fin(a L ) , , ( -hn{aL,b L ) 1 



mi ^VmKMf-^) (m|-m|) 

2f2n(a L ) ^ 1 (hn{a>R) hn(a L ) 



mi (m 2 R -m 2 L ) V m| m| 

where = Mf 2 /m 2 L f° r ^(1) or SU(2) contributions, respectively, while 
or = M\jm\ and = ^ jm 2 LR . The loop functions /j( Ci „)(a;)'s are reported 
in appendix A while fi( c , n )(x,y) = fi( c , n )(x) - /i( c , n )(y) and t /3 = tan/3. 
In the above equations, for completeness, we included the subdominant con- 
tributions that are not tan /3 enhanced. 

In fig.l we can see the bounds on 8 2 ^ and 5g 2 L coming from \x — > &-/ and 



LL 
32 



r — > /i7, respectively. As for r — > e7, the limit on 5^ is such that 5^/5. 
(Br(r — ► e^f)/Br(r — > ytry)) 1 / 2 . The red lines correspond to the full computa- 
tions while the green and the blue lines refer to the mass insertion (MI) and 
to the generalized mass insertion (GMI) approximations, respectively. 
In the case, MI and full computations give indistinguishable results and 
the GMI gives a very satisfactory approximation (see the next section for a 
quantitative estimate). In the 5^ 2 case the degree of approximation is worse 
than in the b 2 \ case. The motivation is that the flavour violating (FV) and 
conserving (FC) insertions of the 32 sector are generally larger than those 
relative to the 21 sector. For instance, 5^ 2 — I0 2 5 2 i and <5|| i = m T /m^5 22 L . 
A large FV insertion, as for 8^ 2 = 0.3, produces a sizable distinction between 
approximated and full computations as it must be when we go away from the 
perturbative region. 

The effect of a large FC insertion is evident for small mo, where the \x term 
is much larger than the slepton masses, so that to treat m T /i tan f3/m 2 LR as 
insertion is not properly correct. 

In fig. 1, the most appreciable deviations between the full computation and 
the GMI one are due to the approximations in the neutralino mixing and they 
tend to vanish in the large tan/3 limit, where the neutralino (chargino) mass 
matrix expansion is better justified. 

We remark that we don't see these deviations in the regions with chargino or 
pure-Bino dominance (in fact these last contributions are lesser affected by 
the GMI approximations). The bounds are rather insensitive to the sign of \i. 

6.2 Bounds on 5lr 

In this case, the only contribution arises from the B exchange and the 
amplitude does not contain a tan/3 factor so b l [ R bound is tan/5 independent, 
unlike all the other bounds. 
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Figure 1. Upper limits on 5[ L from Z, — ► h^y. In the plots we set \i > 0, tan/3 = 10 
and a=0. Red lines correspond to full computation, green and blue lines to MI and 
GMI approximations respectively. 



In the GMI approximation, the amplitude has the following expression: 



A ij = gl A RL ( Ml\ f hnjaR) _ hnM \ 

L1 4tc (m 2 L — m 2 R ) \mi.J V m 2 R m 2 L J 

We note that the chirality flip is realized directly by the mass insertion so 
we can understand the order of the bound, compared to the LL case, S^ R ~ 
(rrii/rh) tan/3^ L , as confirmed numerically. 

While MI and full computations give practically the same result both in 21 
and in 32 sectors, the GMI approximation starts to work very well when 
Mi/2 > 300Ge\^. This result can be better understood by bearing in mind the 
conditions under which the GMI approximation can be applied. A necessary 
condition is that M 12 > m w which is not satisfied when Mi/ 2 < SOOGeV^. 
In the LL case we did not have this problem because of the chargino domi- 
nance. We remind that M 2 ~ 2M\ so we would expect the same behavior when 
^i/2 < VsQGeV but this region is forbidden by the LSP Bino constraint. 
In the 32 sector, both MI and GMI approximations show a sizable deviation 
from the full computation in small regions because of a large 5 Z RL . 
A large 5^ L induces a mass split of order 8 R \ in the third generation so, in 
the mass insertion language, we are neglecting next to leading terms of order 
0((5 R \) 2 ) not so suppressed in the examined case. We note that the same 
argument holds for the 8 Z2 L case. 
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Figure 2. Upper limits on 5 RL from li — > Ij'y. In the plots we set \i > 0, tan /3 = 10 
and a=0. Red lines correspond to full computation, green and blue lines to MI and 
GMI approximations respectively. 

6.3 Bounds on Srr 



The Srr sector requires a bit more of attention because of some cancella- 
tions occurring among the amplitudes in some regions of the parameter space. 
The bounds on Srr become very weak so, it is interesting to check what is 
limit of applicability of the mass insertion approximations 2 . 
The origin of this cancellations is the destructive interference between the 
dominant contributions coming from the B (with chirality flip implemented 
through a FC mass insertion) and BH° exchange. To better understand the 
nature of these cancellations let us derive, in the GMI approximation, the 
amplitude associated to Srr: 



Ml A ij 



4/ln(flfl) 
/ 2/2n(flfl) 



+ 



+ 



III 



R 



(rnj j —m 2 R ) \ m 



{m\-m\ 



In the above equation, the first term is relative to the pure Bino amplitude 
with external chirality flip, the second one corresponds to the BH° mixing in 
the neutralino sector while the last term originates from the B contribution 
with internal sfermion line chirality flip. 

2 This typically occurs for RR type MI as long as universality in the gaugino masses 
is maintained at the high scale. Although in a completely generic situation without 
any universal boundary conditions, such cancellations can occur for LL type MI also 
[36]. 
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It is easy to check numerically that the dominant contributions, proportional 
to tan /3, have opposite sign in all the parameter space. To get a feeling of 
the reason of the above opposite sign, we note that the amplitude relative to 
B (with chirality flip implemented through a FC mass insertion) is propor- 
tional to NllNrr with sign(NLLNfiR) = — 1 while the contribution arising 
from the BH° exchange is proportional to N LR N RR with sign(N LR N RR ) = +1 
(see eqs.3). Vice- versa, the same type of contributions in the Sll case have 
the same sign being proportional to N RR N LL with sign(N RR N LL ) = —1 and 
to N RL N LL with sign(N RL N LL ) = —1, respectively. The above difference de- 
pends on the opposite sign between the hypercharge of SU(2) doublets and 
U(l) singlets. 

If some cancellations occur among the leading contributions, subleading ef- 
fects, generally disregarded, could become important or even dominant. 
In this spirit, we retain the amplitude relative to a chirality flip realized in the 
external fermion line, neglected in [14] because it is not tan (3 enhanced. 
Moreover, the above amplitude shows that while the dominant (tan/5 en- 
hanced) contributions are proportional to the \x mass term 3 the pure B am- 
plitude with external chirality flip is fi independent . In this way, the branching 
ratio is not invariant under the change of the \x sign so, in general, one has to 
consider both cases. 

In fig. 3 we show the upper limits on 5 2 RR from fi — > cy. The limit on 8^ R are 
simply obtained by S^/S^f 1 = (Br(r — > lj^)/Br(fj, — > e'j)) 1 ^ 2 thus, by now 
<5|| R and are not constrained at all. 

As we can see, we are not able to remove these cancellations, in fact, the effect 
of the tan j3 independent contribution is only a shift of the cancellation region 
(the same thing happens if we flip the \i sign). 

We remark that, such cancellations occur in all the situations where the con- 
tribution of the A term to 5 2 RL = (A — //tan^m^ is negligible. There are 
well known model dependent upper bounds on the A parameter to avoid color 
and e.m. charge breaking, in particular \A\/m R < 3 in mSUGRA. Moreover, 
mSUGRA requires a large \x term to fulfill the electroweak symmetry break- 
ing, thus, we cannot invert the relative sign between the two amplitudes. In 
this spirit, we neglected the A term in A\. It is noteworthy that the MI ap- 
proximation works very well reproducing the same cancellation regions as the 
full computation. In the GMI case, we have a net shift of this region but the 
general structure is maintained. 

In conclusion, /i — > doesn't allow to put a bound in the RR sector, so, we 
take into account other LFV processes as fi — > eee and /j, — > e in Nuclei. 
As fig. 5 shows, we find that the last process suffers from a bigger cancellation 
problem than \i — > (in the RR sector) while \i — > eee does not. 
This result requires some explanations. As we have seen in sec. 5, the dom- 
inant contribution to Br(fi e in Nuclei) and Br(fi — > eee) arises from 

3 In reality, this is true only if we neglect the A term in 5 2 R 2 L = (A — //tan /3)m M in 
the pure B amplitude. 
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Figure 3. Upper limits on 5jl R from \i — ► ej; In the plots we set \i > 0, tan/3 = 10 
and a=0 and red lines correspond to the full computation, green and blue lines to 
the MI and the GMI approximations, respectively. 

the dipole operator (that is tan/3 enhanced) but, if the dipole amplitude is 
strongly suppressed by some cancellations, Br{fi — > ej) goes to zero while 
Br {fx — *> e in Nuclei) and Br{fi — > eee) are dominated by non-dipole con- 
tributions. So, in principle, fi — > e in Nuclei and fi — > eee could be able to 
bound 5rr. As we can see in fig. 5, this is the case of fi — > eee that gives a 
bound for Srr < 0.4 that is, correctly, tan/3 independent. 
On the other hand, we find that Br{fi — > e in Nuclei) has additional cancel- 
lations between dipole and not-dipole amplitudes. As a final effect, we have 
that fi — > e in Nuclei suffers from a similar cancellation problem as ft — > e'j 
but, as can be expected, in a different region of the parameter space. 
Because of this strong cancellations, fi — > ej and fi — > e in Nuclei prevent us 
from getting a bound in the RR sector both at the present and even in the 
future when their experimental sensitivity will be improved. 
However, an interesting feature is that u —>■ ej and fi — ► e in Nuclei ampli- 
tudes have cancellations in different regions, so, if we combine the two pro- 
cesses, we obtain a more stringent bound {S^f 1 < 0.2) than the one coming 
from fi — > eee {d^f 1 < 0.4). It is noteworthy that the study of combined pro- 
cesses allows to extract additional informations respect to each separate case. 



6.4 Bounds on the double mass insertion 5 32 5 31 from fi —>■ 

As we have seen in the previous section, the bounds on r — > /ry and 
t — > e7 are very loose. This is due to a worse experimental resolution on 
the above processes compared to the fx — > e^y one. However, we can extract 
additional information in the 32 and 31 sectors applying to fi — > erf. 
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Figure 4. Upper limits on 5^ R from u — > e in Ti (a) and /i — > eee (b). In the plots 
we set fi > 0, tan/3 = 10 and a=0. 




Figure 5. Branching ratios of fj, — * e transitions normalized to the actual experimen- 
tal upper bounds vs S 2 R l R . Red dots correspond to X = 7, blue dots to X = "in Nuclei" 
and green dots to X = ee. In the scatter plots we have taken /i > 0, 3 < tan (3 < 40 
50 GeV< m <500 GeV, 50 GeV < M 1/2 < 500 GeV and -3 < a < 3. On the right, 
for each point of the parameter space, it is shown only the biggest normalized Br 
value among \i — > ej, \i — > eee and fx — > e in Ti. 

The point is that we are able to put bounds on the product of two mass 
insertions, namely 532^31 • In general, we can pass from the second to the first 
generation or through the 821 or through the 523^31 insertions. 
Now, to constraint the Mis, we proceed exactly as for a single MI. We put 
to zero all off diagonal FV entries in the slepton mass matrix except for the 
two Mis we are interested to bound. At this point, we give the analytical 
expressions for all the ^23^31 type insertions in the GMI approach. 
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For S^S^i we get the following amplitude: 
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The amplitude A%\, relative to S^fS^f-, is obtained by A R \ = A 2 L \(L^R). In 
fig. 6 we show the bounds on S^S^ and on S^S^f 1 . 

As we can see, they exhibit different behaviors, especially for m smaller than 
M1/2 due to the m R and m L mass difference (in fact, while in the first case we 
have two left handed and one right handed sfermions running in the loop, in 
the second case we have the opposite situation). 

So, while 5 RL and 5 % [ R are indistinguishable, it is not so for #2*3^31 and S^fS^. 
The amplitude associated to S^S^i, namely A\\ = (A L \)su(2) + (^-1,2)^(1), 
reads: 
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The contribution arising from a S^ L 5 RR type insertion reads: 
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In the fig. 7 we show the bounds for #23 "31 and for #23 "31 



23$3i R (equal to the 
bounds on S^S^). It is to note that S^S^ is strongly constrained because, 
the associate amplitude, is m T /m^ enhanced respect to the usual Bino-like me- 
diated processes (being the chirality flip implemented in the internal sfermion 
line through 8^ and not by 5 22 R , as usual). The amplitude A 2 R \, relative to 
5« R 5^, is A 2 R \ = A? 3 (L~R). 

Finally we derive the expression for the amplitude associated to S 23 R S§^: 
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We note that, in general, #23^31 has a bound comparable to 5 2 i, being 8 2 i — 
^23^31- The only exception is for S^S^ (S^S^f ) as discussed above. 




Figure 6. Upper limits on 5 R \ 6j^ L and S^ R 6 RR from \i — » erf. We have set fi > 0, 
tan/3 = 10 and a=0. Red lines correspond to full computation, green and blue lines 
to MI and GMI approximations respectively. 




Figure 7. Upper limits on 5^ L S RR and 5f? L S^ L from \i — ► e'y. We have set \i > 0, 
tan (5 = 10 and a=0. Red lines correspond to the full computation, green and blue 
lines to the MI and to the GMI approximations respectively. 
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Figure 8. Ratios between exact and MI approximation (red dots) and between exact 
and GMI approximation (blue dots) vs 5 l l L . On the right, tiny light blue dots are 
relative to the chargino-neutralino approximations while tiny blue and red dots are 
due to large FC terms effects. 

7 Mass eigenstate vs mass insertion: a numerical analysis 

At this point, our purpose is to make a quantitative estimate of the good- 
ness of the MI and GMI approximations. 

Even if we have already seen that the above approximations give us the same 
bounds as the full computation, we were not able, in the previous analyses, 
to quantify and to distinguish correctly the approximations induced by the 
chargino (neutralino) branch and by the slepton branch. In the slepton case, 
we can distinguish between two sources of approximations, namely, FC and 
FV terms. 

In fig. 8 we show the ratios between full/Ml computations (red dots) and 
full/GMI computations (blue dots) vs. FV mass insertions. 
As we can see in the 21 sector, MI approximates correctly (at (5—10)% level) 
the full computation until large FV insertions (^i — 0.2), while the two com- 
putations are practically indistinguishable for 5 2 i < 0.1. The induced approx- 
imations are easily understood if we remind that the approach we followed to 
put the bounds on the 5 l i insertions was to consider only one mass insertion 
contributing at a time for the Br(l{ —>■ Ij'y)- If we stop in the first term we 
neglect terms of order 5fj in the amplitudes inducing, naively, an approxima- 
tion on the branching ratios (Br(Sij + 5fj) — Br(5ij))/Br(5ij + 5fj) ~ 25?-, as 
it is well reproduced numerically. The FC insertion 5^ does not produce any 
sizable effects, in fact, in the worse case (for large /1 and tan/3 and moderate 
slepton masses), we have S^p < 0.1. The last argument is not true in the 32 
sector where, being S^/ 5^ — mr/rrin, we can have a not perturbative FC 
insertion. This is clear in fig. 8 where, the MI and the GMI approaches under- 
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line a (40 — 50)% deviations respect to the exact calculation (tiny dots refer 
to the tan/? > 30 and My 2 > 300, region where we have sizable deviations 
due to the slepton FC terms only). 

Now, we want to discuss the approximations brought from the chargino and 
the neutralino sectors in the GMI approach (the following discussion is obvi- 
ously flavour independent). 

As discussed in section 3, it is allowed to use this method when the elements 
outside from the diagonal (proportional to m w ) are much smaller than those 
diagonals (// and M 2) i). 

On the left of fig. 8, tiny blue dots show a (20 — 30)% approximation (it hap- 
pens in M\/2 < 300 and tan/3 < 10) while the much larger ones refer to the 
M\/2 > 300 and tan/3 > 10 region where the GMI conditions are fulfilled. 
In the 32 case light blue dots correspond to deviations induced by the GMI 
approximation in fact they are related to a range of parameters (m > 300 
and tan/3 < 5) where 8^ < 0.1. 

To summarize the results found, we can say that MI approximation produce 
the same features as the exact calculation even if strong cancellations occur. 
The approach works better in the 21 sector than in the 32 one because, while 
in the first case we always have perturbative FC terms, in the second case it 
is not so and we can reach sizable deviations (up to a (40 — 50)% level) from 
the full computation. 

Moreover, we have a 10% approximation until FV terms of order 0.2. 
The GMI approximation works very well, as the MI approximation, up to 
gaugino masses heavier than l^OGeV and it produces the cancellations in a 
shifted region respect to the exact case. 

In conclusion we can say that, except for fine-tuned cases, the last approach 
is very satisfactory. 



8 Conclusions 

In this work, in the first stage, we studied the constraints on flavour vio- 
lating terms in low energy SUSY coming from Zj — > ^7. 

We have carried out the analysis both in the mass eigenstate and in the mass 
insertion approximations clarifying the limit of applicability of these approxi- 
mations. 

In particular, we focused on the RR sector where strong cancellations make 
the sector unconstrained. 

We showed that these cancellations prevent us from getting a bound in the 
RR sector both at the present and even in the future when the experimental 
sensitivity on Br(k — > ^7) will be improved. 

Finally we took into consideration the bounds on the various double mass 
insertions: 8^8^, 8^8^, 8^8§[ R , 8^8^ and 8^8^. It is clear that the 
bounds are approximatively the same as in 8 2 \ being 8 2 i — ^23^31 except for 
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the last one suppressed by a m^/m T factor. So, in spite of very weak bounds 
on and on # 31 coming from Br(r — > /17) and Br{r e^y) respectively, we 
have stronger bounds on their product thanks to Br{ji — > ej) experimental 
sensitivity. These limits are important in order to get the largest amount of 
information on SUSY flavour symmetry breaking. 

Summarizing the results found in this first stage, we can say that MI ap- 
proximation produce the same features as the exact calculation even if strong 
cancellations occur. 

The approach works better in the 21 sector, up to (5 — 10)% approximation 
level until 8 2 i < 0.2, than in the 32 sector, where, large off diagonal flavour 
conserving terms can induce a (40—50)% deviation from the full computation. 
The GMI approximation works very well, as the MI approximation, except for 
some special regions where induces a (20—30)% approximation with respect to 
the exact calculation and produces the cancellations in a shift region respect 
to the exact case. In conclusion, we can say that, except for fine-tuned cases, 
the last approach is very satisfactory. 

In a second stage, being our aim to find constraints in the RR sector, we ex- 
amined other LFV processes as /x — > eee and /x — > e in Nuclei. 
We found that the last process suffers from a bigger cancellation problem than 
/i — > e7 (in the RR sector) while /x — > eee does not. 

However, an interesting feature is that /x — > e^ and /x — > e in Nuclei ampli- 
tudes have cancellations in different regions, so, if we combine the two pro- 
cesses, we obtain a more stringent bound (5^ < 0.2) than the one coming from 
/x — > eee (S^ 1 < 0.4). It is noteworthy that the study of combined processes 
allows to extract additional information respect to an individual analysis of 
all these processes. In particular, it makes it possible to put bounds on sectors 
previously unconstrained by /x — > ej. 
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A Loop functions 



In this appendix we report the explicit expressions for the loop functions 
appearing in the text: 

-17a 3 + 9a 2 + 9a - 1 + 6a 2 (a + 3) In a 
hn[a) ~ 24(1 - af 

-5a 2 + 4a + 1 + 2a(a + 2) In a 
hn[a) = 4(1 - a)* 
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fic(a) = 
Jin (a) = 



1 + 2a In a — a 2 
hn[a) = 2(1 -a) 3 
—a 3 — 9a 2 + 9a + 1 + 6a(a + 1) In a 
6(1 - a) 5 
-a 2 - 4a + 5 + 2(2a + 1) In a 
= 2(1 - a) 4 

3a 4 + 44a 3 - 36a 2 - 12a + 1 - 12a 2 (2a + 3) In a 



hc{o) 



24(1 - a) 6 
a 3 + 9a 2 — 9a — 1 — 6a(a + 1) In a 



4(1 -a) 5 

10a 3 + 9a 2 - 18a - 1 - 3a(3 + 6a + a 2 ) In a 



hr (a) 



6(1 -a) 6 
3a 2 - 3 - (a 2 + 4a + 1) In a 
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